Spin and the Honeycomb Lattice: Lessons from Graphene 
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Spin-1/2 particles such as the electron are described by the Dirac equation, which allows for two 
spin eigenvalues (up or down) and two types of energy eigenvalues (positive or negative, correspond- 
ing to the electron and the positron). A model of electrons hopping from atom to atom in graphene's 
honeycomb lattice gives low-energy electronic excitations that obey a relation formally identical to a 
2+1 dimensional Dirac equation. Graphene's spin equivalent, "pseudospin" , arises from the degen- 
eracy introduced by the honeycomb lattice's two inequivalent atomic sites per unit cell. Previously 
it has been thought that the usual electron spin and the pseudospin indexing the graphene sublattice 
state are merely analogues. Here we show that the pseudospin is also a real angular momentum. 
This identification explains the suppression of electron backscattering in carbon nanotubes and the 
angular dependence of light absorption by graphene. Furthermore, it demonstrates that half-integer 
spin like that carried by the quarks and leptons can derive from hidden substructure, not of the 
particles themselves, but rather of the space in which these particles live. 
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"Spin" refers to an angular momentum that has no 
classical analogue — it is not possible to understand spin 
in terms of a mechanical model of some rotating object 
The net angular momentum of composite particles, 
such as protons, neutrons, atoms, and molecules, derives 
from the spins of their constituents, plus any orbital an- 
gular momenta due to the constituents' relative motion. 
Spin and orbital angular momenta are quantitatively dis- 
tinguishable, since the former can be half-integer while 
the latter take on integer values only, measured in units 
of the reduced Planck constant h. In the standard model 
of particle physics the ultimate constituents of matter are 
the quarks and leptons. These particles have no internal 
structure down to length scales of 10 -18 m (limited by 
the collision energies ~ 200 GeV currently achievable in 
particle accelerators) 0, so their spins are considered 
intrinsic. 



The deepest insight into the origin of spin has been 
provided by Dirac, who manipulated Einstein's quadratic 
energy-momentum relation E 2 = p 2 c 2 + m 2 c 4 to give a 
linear equation consistent with the postulates of quantum 
mechanics Dirac's equation predicts not only spin 
but also antiparticles, which were unknown at the time. 
Thus we understand, for example, the electron's spin h/2 
and the existence of the positron as natural consequences 
of the Dirac equation, which is built on the theories of 
relativity and quantum mechanics. As was pointed out 
25 years ago [H, H|, the low-energy electronic excitations 
in graphene obey a 2 + 1 dimensional Dirac equation, 
with holes and the sublattice state playing the role of 
positrons and spin respectively. In this Letter we show 
that the sublattice state vector describes an 'intrinsic' an- 
gular momentum in 3 + 1 dimensions. This identification 
provides a physical model that associates spin with an 
underlying structure. Unlike the case of composite parti- 



cles, where spin follows from other spins, in this example 
the spin h/2 is a consequence of the nontrivial spatial 
lattice, invisible at low energies, that hosts the particle. 

Graphene's Dirac equation follows from the tight- 
binding (i.e. hopping) model, which was first applied to 
graphene by Wallace M ■ With the experimental isolation 
of carbon nanotubes (71 and, more recently, graphene it- 
self Q, the hopping model has been shown togive an 
effective description of these real materials 0, [Iq| . Fur- 
thermore, graphene equivalents of the quantum relativis- 
tic effects implied by the Dirac equation, such as Klein 
tunneling and Zitterbewegung, are now experimentally ac- 
cessible, making this condensed matter system a practical 
test bed for these particle physics phenomena [ill Il2| • 

We first show how the hopping model produces the 
2+1 dimensional Dirac equation, working in 3+1 dimen- 
sions but without initial reference to specific coordinate 
axes. Graphene's electronic states are described as a lin- 
ear combination of atomic orbitals constructed to fulfill 
the Bloch condition *Q(r + R) = e iQ ' R *Q(r), 



*Q(r) 



3 



N 



Q' 



(r-R/) + c^ 



(1) 

where j — (m, n) indexes the N sites of the hexagonal 
Bravais lattice described by Rj = mai + n&2, and the 
vectors Rj 1 (Rf ) point to the 'A' ('B') sublattice sites 
within the unit cell j, respectively (see Fig. [T]). Indices 
labeling the usual electron spin have been suppressed for 
notational convenience and will be henceforth. The co- 
efficients cq multiplying the carbon atoms' 2P Z atomic 
orbitals (f>(r) are chosen to solve the crystal Hamiltonian. 
In the literature there are equivalent alternatives [H, @, H3] 
to the expansion ([I]), but this choice reveals the momen- 
tum space symmetry of the Hamiltonian more readily 
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FIG. 1: The honeycomb lattice in real space (left) and the 
corresponding band structure in reciprocal space (right). On 
the left, gray and black circles represent A and B lattice sites 
respectively. One choice of unit cell is demarcated with a 
dashed line. Coordinate axes oriented to give the Pauli and 
Dirac forms of the electronic Hamiltonian are labeled "P" and 
"D" . On the right, the hexagonal first Brillouin zone is shown, 
with the positions of K + points indicated by thin arrows. 
While generally the contours of constant energy [based on 
the full Hamiltonian Q] show threefold or sixfold rotational 
symmetry, near the points they are circular [a. fl3|. 



The electronic Hamiltonian contains two kinds of 
terms, one representing an electron's energy on a par- 
ticular site, and the other representing the energy ad- 
vantage conferred by the freedom to hop to a neighbor- 
ing site. The amplitude for nearest-neighbor hopping is 
parametrized by an energy t, and next-nearest-neighbor 
and higher order effects are neglected. An electron occu- 
pying the state (f>{r— Hf) is described by an operator A^ R . 
that creates an electron on the "A" site in cell j when it 
acts on the vacuum state |0). With similar language for 
the "B" sites, the graphene Hamiltonian is 



H = 



j 



<i,3> 



(2) 

where the notation indicates sums over all sites j 

and their nearest neighbors i. In graphene the A sites 
are identical to the B sites modulo a ir rotation, so the 
energies £a = £b- Following Semenoff we allow for 
different site energies (£a ^ £b), as in hexagonal boron 
nitride, and take the graphene limit where appropriate. 

Introducing the Fourier transform of the operators 
^R 4 = V , Aq expi/'R, • Qj)/y/N, where the Qj = 
Tjr-bi + jj^i>2 are the N = N1N2 wave vectors in the 
first Brillouin zone, puts the Hamiltonian ([2]) in the form 



H 
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H 



(3) 



where we have used the closure relation ^T. e tRj ' ^ ^ 
NSq qi, and defined the single particle Hamiltonian 

-A/t 1 + e - iQai + e - iQa2 

1 + e jQai + e iQa2 A/t 



n = -t 



Here the energy difference A is defined by A = (£a — 
£b)/2, and the energy origin is chosen such that (£a + 
£b)/2 = 0. Graphene has two atoms per unit cell, each of 
which donates one electron to the valence band, so in the 
lowest approximation the first Brillouin zone is exactly 
filled. Thus the Fermi energy is zero here. 

The off-diagonal matrix elements of the T~L matrix (|4]) 
vanish at the corners Q = K of the Brillouin zone (see 
Fig. QJ, which in the graphene case £a — £b leads to 
the famous degeneracy at these 'Dirac points': K K = 
mbi + n\>2 13]. Here n = ±1 is the "valley" 
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index that distinguishes the two inequivalent types of K 
points. To see the structure of low-energy excitations 
near the Dirac points we define k = Q — K and restrict 
our analysis to the case where k • is small. Recalling 
a; • hj = 2-jrSij , to lowest order 



H 



A 

\/?>tak.{Hihd 



\Z3tak(K&d - 
ik s )/2 -A 



ia s )/2 



(5) 



where we have defined difference and sum unit vectors 
a^ = (ai — a 2 )/a and a s = (ai + a 2 )/V3a. (A third 
vector a„ = a<j x a s normal to the plane will also prove 
useful.) The Hamiltonian (J5j> is a rotational invariant, 
depending only on scalars and the scalar products of 3D 
vectors. As shown in Fig. [TJ the eigenvalues of H are 
independent of the direction of k near K, which is to say 
that the Hamiltonian is effectively isotropic in the plane 
— the lattice has disappeared. 

The analysis leading to ([5]) has been entirely in terms 
of the basis vectors of the direct and reciprocal lattices 
of the honeycomb structure, and is independent of any 
choice of orientation for a 3D Cartesian coordinate sys- 
tem. The predictions of the theory are independent of 
this choice. Two particular orientations put % in fa- 
miliar and suggestive forms. Choosing the coordinate 
orientations labeled "D" and "P" in Fig. [TJ defining the 
Fermi velocity vf = V3at/2h, and writing the conjugate 
momentum p = fik gives 

T-Ld = v F {na x p y - a y p x ) + a z A and (6a) 
Hp = v F (na x p x + cFyPy) + a z A, (6b) 

where the aj are the usual Pauli matrices. The matrix 



(4) 



% d (|6a|) illustrates the connection between the graphene 
Hamiltonian and the Dirac equation. In 3 + 1 dimen- 
sions the Dirac Hamiltonian is H = C7 (7 • p + mc), 
with 4x4 matrices 7 M = (7°,7). In 2 + 1 dimen- 
sions the necessary anticommutation relations {7^, 7"} = 
2g tJ,u = 2 x diag(l, — 1, — 1) can be satisfied with a 2 x 2 
representation such as 7 M = (7 ,7) = (a Zl ia x ,iKa y ). 
With this definition the Hamiltonian matrix (|6a|) be- 
comes Ti. = fF7°(7 • P + rn'vF), with an effective mass 
m! defined by A = m'vp. (We designate the first 
two components of a 3D vector with an arrow, e.g. 
p = p x x + p y y + p z z = p + p z z.) Thus this hopping 
model gives low-energy electronic excitations that obey 
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a 2 + 1 dimensional Dirac equation, with the Fermi ve- 
locity vf playing the usual role of the speed of light c. 

With the P coordinate orientation graphene's Hamilto- 
nian matrix has the convenient abbreviation % = vpa -p 
near K + . Sometimes this is written % — vf<j ■ p, which 
can give the impression that either a z or p z is strictly 
zero. In fact, both cr and p have three nonzero compo- 
nents; it just happens that a z andp 2 do not appear in the 
graphene W. While any state confined to the sheet must 
have an expectation value (p z ) = 0, because the electrons 
are localized in z the characteristic magnitude of p z is 
h/a z , where a z is a length scale measuring the z extent 
of the 2P Z orbitals of the expansion ((T|) . Taking p z = is 
just as improper for an electron in a honeycomb lattice 
as it is for an electron in an atomic Coulomb potential. 
The case of a z will be taken up below. 

It is not obvious whether the operator cr, which indexes 
the "pseudospin" arising from the AB sublattice degen- 
eracy, corresponds to a real angular momentum 0, Il4j . 
It might describe a merely analogous two-state system, 
borrowing the same SU(2) algebra like the isospin sym- 
metry connecting the proton and neutron. Comparison 
with the 3+1 dimensional Dirac equation makes this sec- 
ond option look likely. The 4x4 Dirac matrices give each 
state two qubits, one each for particle or antiparticle and 
spin-up or spin-down values. Since the pseudospin labels 
the band index /? [see Eq.|H] below] and the 2D Dirac ma- 
trices are only 2 x 2, it would seem that the one qubit 
available is engaged. 

Surprisingly, this one qubit manages to cover both vari- 
ables. The pseudospin is related to a real angular mo- 
mentum, as can be discovered by calculating the com- 
mutator of the Hamiltonian with the orbital angular mo- 
mentum L = rxp (Ref. Q), 

/ VyPz \ 

[H P ,L] = -ihv F \ -na x p z . (7) 

\K(T x Py - (TyPxJ 

That the Hamiltonian has rotational symmetry about 
the axis perpendicular to the plane and the commuta- 
tor [Hp,L z ] is not zero together indicate that there is 
another angular momentum in the problem. 

In coordinate-independent notation the honeycomb 
Hamiltonian ([5]) is 

% = • u, where (8a) 

h 

S = ^(kct x &d + o y a s + kct z a„), and (8b) 

u = (p ■ k d ) a d + (p • a s ) a s + k(A/v f ) a„. (8c) 

The operator S defined by (|8b|) we term "lattice spin" 
to distinguish it from both the dimensionless pseudospin 
and the usual electron spin. One can easily verify that 
[H, (L + S) • a n ] = for any value of A. Thus neither 
the lattice spin nor the orbital angular momentum is sep- 
arately a constant of the motion, but the projection of 



D H j[H,r] j[H,L] |[g 1 |oj conserved system 

3 cer ■ p cer cerxp — cerxp J neutrino 

2 vS ■ p vS iiirxp —vcrxp J z graphene 

1 vcfzPz vgz vtJzixp —vcrxpzi L z , Sz nanotube 

TABLE I: Commutation relations for massless Dirac Hamil- 
tonians H. In the Heisenberg representation an operator A 
without explicit time dependence obeys the equation of mo- 
tion dA/dt = j?[H, A]. The system shorthand refers to right- 
handed massless neutrinos, graphene in the xy plane, and a 
metallic carbon nanotube with its axis along the z direction 
respectively. 



the total angular momentum J = L + S onto the plane- 
normal axis (the z-axis for P or D coordinates) is a con- 
served quantity. It is not possible to confuse S with an 
orbital angular momentum, as its eigenvalues have half- 
integer magnitude in units of h. 

In comparison, the connection between pseudospin and 
angular momentum in a nanotube is less obvious, for in 
that system the component of the orbital angular mo- 
mentum along the nanotube axis commutes with the 
Hamiltonian. In a carbon (A = 0) nanotube the com- 
ponent of S along the axis also commutes with H, so the 
axial projections of the lattice spin and the orbital angu- 
lar momentum are conserved separately. Table U summa- 
rizes the relevant angular momentum commutation rela- 
tions for Hamiltonians of the type H — u X^f a iPi with 
number of spatial dimensions D = 1, 2, or 3. 

Since earlier work asserts that the pseudospin is not as- 
sociated with an angular momentum 0, [3] , it is worth 
exploring why we are led to a different result. In a strictly 
2D system angular momentum is defined only in a lim- 
ited sense. There is only one generator corresponding to 
(commuting) rotations about the direction normal to the 
plane [13], which is inconsistent with a 3+1 dimensional 
spin. By confining the electrons, the graphene sheet pro- 
duces an electronic Hamiltonian % oc S • u which lacks 
the full 3D rotational symmetry of the vacuum. However, 
the (scalar) Hamiltonian is invariant under rotations of 
the system relative to the three coordinate axes, and the 
generators of these rotations, i.e. the angular momentum 
operators, are still well defined. (Table [TTI juxtaposes the 
rotational properties of various spin-type Hamiltonians.) 
And unlike two-level systems that can be reduced to the 
spin- 1/2 problem by analogy only, the honeycomb struc- 
ture produces a spin variable specified in relation to a 
laboratory coordinate system. In other words, the direc- 
tion of u, like that of a magnetic field B, can be related to 
'east', where no such relation appears in, say, the isospin 
problem. Since the operator i?(n, $) = exp(— i$S ■ h/h) 
generates a rotation of the observable S by an angle $ 
about the n = n/|n| axis, where n is a direction in real 
space, S must be a real angular momentum. See the 
Appendixes for a proof and further discussion. 

The analogy between the vector u and a magnetic field 
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neutrino [15] 


spin-orbit [16] Z 


eeman[16l honeycomb deuteronfl. 171 


Hamiltonian H oc S ■ Z 


Sp 


S • L 


S B 


S ■ u 


S (i). S W 


S name 


spin 


spin 


spin 


lattice spin 


isospin 


3D rotational symmetry 


yes 


yes 


no 


no 


yes 


[H,J] oc 








S x B 


Ci 





[n • J,S]/ih 


Sxn 


Sxn 


Sxn 


Sxn 





[n • J, Z]/iH 


pxn 


Lxn 










S = angular momentum 


yes 


yes 


yes 


yes 


no 



TABLE II: Hamiltonians with spin-like variables S. A general unit vector is designated by n, and the placeholders are 
£i = S x (u — p) + (S • a„)(a„ x p) and $2=px While all Hamiltonians are rotational invariants in the 

sense that they are scalars, the Zeeman and honeycomb Hamiltonians do not represent rotationally invariant systems. Also 
note that B, u, and the isospin S are not spatial, quantum mechanical vectors V, as they do not satisfy [n-J, V] = iKV x n.[l6l] 
However, as the Zeeman example illustrates, this condition is not required for S to be an angular momentum. 



B is helpful for understanding the close relationship be- 
tween energy and angular momentum that follows from 
the er operator's coverage of both qubits. We compare 
the time evolution operator U = exp(— iHt/K) with the 
rotation operator R(u, $). For the honeycomb Hamil- 
tonian (|8a|) these expressions are identical; time evolving 
the state is equivalent to rotating it about the u axis at a 
rate d<&/dt = (E+ — E-)/h. The energies Ep correspond 
to the two eigenvalues of the Hamiltonian H, 



E 



(9) 



where we have defined the band index ft = ±1. An equiv- 
alent relationship between time evolution and rotations 
appears when a magnetic dipole is placed in a magnetic 
field B, where it is called Larmor precession. The time 
evolution of the lattice spin operator (see TableU]) is given 
by 



dS 

~dt 



2v F 



uxS, 



(10) 



which is exactly the equation of motion of a mag- 
netic dipole in a magnetic field, with the substitution 
2vfu/H — > — 7B (here 7 is the gyromagnetic ratio). 
Thus the mixing between helicity eigenstates produced 
by a mass term in the 2+1 dimensional Dirac Hamilto- 
nian is formally identical to the mixing between spin-up 
and spin-down states produced by perpendicular mag- 
netic field Bj^. In the general case of u not parallel to 
S, Eq. (fTUj) shows that the components of S perpendic- 
ular and parallel to the plane interconvert as a function 
of time. 

Relating pseudospin to angular momentum provides 
an intuitive explanation, the existence of which has long 
been suspected [l8| , for some properties of graphene and 
related materials. For instance, in metallic nanotubes 
the backscattering of electrons is suppressed if the scat- 
tering potential has a range that is long compared to a 
lattice constant [HI, E3] • Since such a potential will treat 
A and B sites identically, it is a scalar with respect to 
the lattice spin, and thus cannot create the lattice spin 



flip required to reverse the electron motion. Conserva- 
tion of pseudospin or lattice spin also provides a way to 
understand Klein tunneling in graphene 12[ , and empha- 
sizes the equivalence between the condensed matter and 
particle physics systems. 

As another example, interpreting the pseudospin as 
connected to a real angular momentum gives a satisfying 
picture of photon-mediated electron-hole pair production 
(or recombination) in graphene. In this process an elec- 
tron with pseudospin parallel to its momentum transi- 
tions between the negative energy valence band and the 
positive energy conduction band, absorbing (or emitting) 
a photon and flipping its pseudospin [20ll2lj. Both bands 
are derived from atomic 2P Z orbitals, and the electron 
spin does not flip, so the usual sources do not contribute 
to the photon's angular momentum. A full description 
[22] ] of this process based on the Hamiltonian (I6bp and 
the substitution p — > p—qA/c (with the quantized vector 
potential A describing the photon) shows that the pho- 
ton's polarization couples the initial and final pseudospin 
states. As the photon's vector polarization is associated 
[lH with its spin h, the form of the transition matrix ele- 
ment requires that the pseudospin be associated with an 
angular momentum h/2. 

Thus all three components of the graphene pseudospin 
have straightforward connections with angular momen- 
tum: reversing an in-plane component creates a spin-1 
particle, while the out-of-plane component clearly gen- 
erates rotations. These complementary examples and 
the interconversion described by the equation of motion 
(|10p make a compelling intuitive case for the pseudospin- 
angular momentum identity. 

Spin states have been previously associated with spe- 
cific types of sites in a lattice within the "stag gere d" for- 
mulation of lattice quantum chromodynamics [231 ] . How- 
ever, the hypcrcubic lattice chosen there is a convenient 
framework for discretizing the Dirac equation, not the 
basis of a more fundamental model from which the Dirac 
equation emerges naturally. The hypercubic lattice can- 
not be viewed as a discretized space through which par- 
ticles move, since necessary phase tuning in the lattice 
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Hamiltonian 23J, |24| ruins its hopping interpretation for 
D > 1. Furthermore, the number of sublattices must be 
fixed artificially using foreknowledge of the Dirac equa- 
tion. These differences are not solely aesthetic; unlike 
the unphysical Hamiltonian contrived for the hypercubic 
lattice, the honeycomb Hamiltonian can be studied in op- 
tical lattices 25[ and natural materials such as graphene 

By generating an internal quantum number that re- 
flects an energy scale much larger than the host particle's 
mass [2] . the graphene example invites connection to the 
intrinsic spin carried by the quarks and leptons. How- 
ever, it is not obvious that a lattice can be found that 
naturally produces the 3 + 1 dimensional Dirac equation 
[26l [27| , let alone the chiral properties of the full stan- 
dard model [28l . |29j . Thus there are two possibilities, 
depending on whether these spins have related origins. 
Either intrinsic spin is also the low-energy signature of 
nontrivial quantized space, or lattice spin represents a 
second, experimentally-accessible type of quantum me- 
chanical angular momentum with no classical analogue. 



Proof 

By definition an operator S is a 3+1 dimensional an- 
gular momentum if its Cartesian components Si obey the 
commutation relations 



[Si, Sj] — iheijkSk 



(11) 



and it generates rotations in real space [16| . The Pauli 
matrices obey [cr^Uj] = 2ieijk<Jk, so S as defined by (|8bj) 
satisfies (fTT|) . The form (|8bj) also makes it clear, through 
the dependence on the a.j, that S is referenced to direc- 
tions in real space. To prove that the Si generate rota- 
tions we consider the Pauli matrices as components of S, 
and rotate S by an angle $ about an axis n according to 
the classical relation |30j 



S' = Scos$ + n(n- S)(l -cos$) + (n x S)sin$, (12) 

and compare the result with the S' found by performing 
the quantum rotation operation 16], 



g/ ^HJ A/h g g — i$J n/?i 



(13) 



For total angular momentum of the form J = S + 
Smother, where [S,J ot hor] = 0, these two transforma- 
tions give the same result, which completes the proof. 



Comparison with the Lorentz group in 2+1 D 

The honeycomb pscudospin er can be associated with 
two distinct algebras: 



[Si,Sj] =ihe ijk S k , {{i,j,k} £ {1,2,3}), 

W,7 v }=2g^, 1/} €{0,1,2}), 

with respective representations 

h. 



r 



^(Ka x ,a y ,Ka z ), and 
(cr z ,ia X: iKa y ). 



and 
(14a) 

(14b) 



(15a) 
(15b) 



As discussed previously, the Si of (I15ap generate rotations 
in 3 spatial dimensions by virtue of satisfying (|14a[) . Sim- 
ilarly, the Y of (fTFb|) are associated with two boosts and 
one rotation in 2 + 1 dimensions by virtue of satisfying 

(EE]). 

The Dirac algebra (|14b[) connects the gamma matrices 
7 M to the Minkowski metric , which here has signature 
(+ — — ). In this representation [3l[ the generators L^ u 
of Lorentz transformations are given by U 1 " = ^[V, r f\. 
Explicit calculation gives 



Ki 


= L 01 = 


— 1 


(16a) 




= L 02 = 


IK 

yo-z, and 


(16b) 


h 


= L 12 = 


K 

2**> 


(16c) 



where, as we will see, a Kj generates a boost along the 
j axis and J3 generates the rotation mixing the {1,2} 
coordinates. 

The form of the Dirac equation must be frame indepen- 
dent. A spinor \x) transforms according to \\') = Dj\x), 
with Dj = exp(—i6Xj) and the transformation generator 
Xj e {Ki, K2, Js}- If a spacetime three-vector, e.g. p M , 
has the corresponding transformation p' M = K^^i^p", 
consistency then requires ^„{j)"f u = Dj 1 -f ,1 Dj. Using 
this relation we find the K^ v (j) : 



/coshfl sinh(9 0^ 
A(l) = sinhfl cosh 9 ) , A(2) = 
V 



^cosh6> sinh# N 

1 
, sinh cosh 1 
(17a) 



'10 
..lid A(3) = ( cos6» -sin9 
, sin 6 cos 9 



(17b) 



which are the two boosts and one rotation expected in 
2 + 1 dimensions. The appearance of hyperbolic functions 
in the Lorentz boosts (I17aj) reveals the distinguishing role 
of the i's which occur in the 7 M representation (| 1 5b[) and 
are notably absent from the Si representation (|15a[) . 
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Note also that while it is possible to define represen- 
tations satisfying (Bl) that do not include k's, the forms 
(fT5)) are preferred. If instead all of the k's were included 
in the definition of u, the expectation value of S for a 
Hamiltonian eigenstate would show undesirable behav- 
ior under change of k. More importantly, including the 
valley index k in the definitions of the Si and the 7^ 
gives these operators the expected transformation prop- 
erties under time reversal, which is implemented by tak- 
ing i — y — i and k — > — K. This prescription can be de- 
duced by noting that % is invariant under time reversal 
and comparing Eq. ((4]) and Eq. ([5]) of the main text. 

Thus the pseudospin Pauli matrices <x can be con- 
nected both to the generators of the Lorentz group in 
2 + 1 dimensions and the generators of the rotation group 
in 3 + 1 dimensions. The choice of algebra corresponds 
to considering the honeycomb lattice as a strictly two di- 
mensional structure, or as a quasi- two dimensional struc- 
ture embedded in three dimensional space. Since ex- 
periments on graphene occur in three dimensional space 
containing inherently three dimensional objects {e.g. the 
photon [22|), this second perspective is sometimes un- 
avoidable. 
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